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A completely integrable nonlinear partial differential equation (PDE) can be associated with 
a system of linear PDEs in an auxiliary function whose compatibility requires that the original 
PDE is satisfied. This associated system is called a Lax pair. Two equivalent representations 
are presented. The first uses a pair of differential operators which leads to a higher order 
linear system for the auxiliary function. The second uses a pair of matrices which leads to a 
first-order linear system. In this paper we present a method, which is easily implemented in 
Maple or Mathematica, to compute an operator Lax pair for a set of PDEs. 

In the operator representation, the determining equations for the Lax pair split into a set of 
kinematic constraints which are independent of the original equation and a set of dynamical 
equations which do depend on it. The kinematic constraints can be solved generically. We 
assume that the operators have a scaling symmetry. The dynamical equations are then reduced 
to a set of nonlinear algebraic equations. 

This approach is illustrated with well-known examples from soliton theory. In particular, it 
is applied to a three parameter class of fifth-order KdV-like evolution equations which includes 
the Lax fifth-order KdV, Sawada-Kotera and Kaup-Kuperschmidt equations. A second Lax 
pair was found for the Sawada— Kotera equation. 
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1. Introduction 

The complete integrability of a system of partial differential equations has been 
a topic of active research over the last forty or so years. Even today there is no 
one generally accepted definition of integrability. Many approaches have been ad- 
vocated (see [32] for a recent review) , all of which have their merits but none seem 
to encapsulate the essence of integrability. However one concept has appeared in 
many different approaches; the Lax pair. This is a reformulation of the original 
system of nonlinear equations as the compatibility condition for a system of linear 
equations. 

The story of Lax pairs begins with the discovery by Lax [31] of such a linear 
system for the ubiquitous KdV equation. This equation [6, 28] models a variety 
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of nonlinear wave phenomena, including shallow water waves [20] and ion-acoustic 
waves in plasmas [1, 4, 11]. Lax's discovery shed light on the then newly discovered 
inverse scattering transform method [16] to solve the KdV equation. It allowed 
this method to be extended to a variety of integrable equations. In 1979 it was 
realised that the Lax pair could be interpreted as a zero curvature condition on 
an appropriate connection [44]. This led to generalisations of the method of in- 
verse scattering transforms. In the 1980s, the algebraic structure of Lax pairs was 
elucidated. The connection between Lax pairs and Kac-Moody algebras was given 
in [13]. Subsequent applications for Lax pairs include Backlund-Darboux transfor- 
mations [35], recursion operators [19] and generating integrable hierarchies via the 
root method [17, 25, 34]. However there was a problem; finding the Lax pair in the 
first place. 

A Lax pair is associated with an infinite hierarchy of local conservation laws; a 
harbinger of complete integrability (though not all equations with a Lax pair have 
an infinite set of conservation laws). Exploiting this connection, Wahlquist and 
Estabrook [14, 40] proposed a method based on pseudopotentials that, in certain 
circumstances, leads to Lax pairs. However their method leads to a nontrivial 
problem of finding a representation of a Lie algebra when only a subset of the 
commutation relations are known. 

Another approach to the construction of a Lax pair is provided by singularity 
analysis. In 1977 it was noted [3] that all symmetry reductions of the classical com- 
pletely integrable equations result in equations that can be transformed to Painleve 
equations. This observation subsequently gave rise to the so-called Painleve test 
[41]. In this approach Lax pairs are generated from truncated Painleve expansions 
[33] . However the complexity of Painleve expansion depends critically on the choice 
of expansion variable. 

Computer algebra packages have also been employed to symbolically verify Lax 
pairs [24, 30] . Of course this requires prior knowledge (or a good guess) of the Lax 
pair. 

In this paper we address the issue of finding Lax pairs, in operator form, for a 
given system of differential equations by an approach that is amenable to computer 
algebra. The determining equations for the Lax pair is split into two sets. The first 
set, the kinematic constraints, are solved generically. The second set, the dynamical 
equations, depend on the system under consideration. We make the assumption 
that the system has a scaling symmetry and that this symmetry is inherited by the 
Lax pair. This allows us to reduce the dynamical equations to an overdetermined 
system of algebraic equations (nonlinear, naturally). Such systems may then be 
solved by Grobner basis techniques. 

In the next section. Lax pairs are introduced in their operator representation. In 
Section 3 the matrix representation of Lax pairs is discussed. The relationship be- 
tween these two representations is also examined. Section 4 introduces the concept 
of scaling symmetries. The algorithm to compute Lax pairs based upon scaling 
symmetries is outlined in the next section. Examples of the algorithm applied to 
some well-known equations are given in Section 6. Finally the algorithm is used 
to classify the integrable subcases of a fifth-order KdV-like equation with three 
parameters. A second Lax pair was found for the Sawada-Kotera equation. 
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2. Leix Pairs in Operator Form 

In this paper we consider nonlinear systems of evolution equations in (1 + 1) di- 
mensions, 

Ut = F(u,Ua;,U2a;,...), (1) 

where x and t are the space and time variables, respectively. The vector u(a;, t) has 
N components Ui and F is a nonlinear function of its arguments. In the examples 

we denote the components of u by u,v,w, Throughout the paper we use the 

subscript notation for partial derivatives. If parameters are present in (1), they will 
be denoted by lower-case Greek letters. 

In his seminal paper [31], Lax showed that completely integrable nonlinear PDEs 
have an associated system of linear PDEs in an auxiliary function ijj{x,t), 

Ctb = Xtb, 

(2) 

il^t = Mi) 

where C and M are linear differential operators (expressed in powers of the total 
derivative operator D-r for the space variable x). il) is an eigenfunction of C corre- 
sponding to eigenvalue A. The operators {C, M.) are now known as a Lax pair for 
(1). The property that (1) is completely integrable is reflected in the fact that the 
eigenvalues do not change with time which makes the problem isospectral. 

Let G be a differential function (functional); that is, a function of x,t,u and 
partial derivatives of u and let Up^q = d^df u. The total derivatives of G is given by 

^ ^ dG ^ dG 

p,<i ' , , 

3) 

dG ^ dG 

p,q P'1 

The sums are finite since we will assume that G depends only on finitely many 
derivatives. 

Example 2.1 The KdV equation [1] for u{x,t) can be recast in dimensionless 
variables as 

Ut + aUUx -I- U3a; = 0. (4) 

The parameter a can be scaled to any real number. Commonly used values are 
a = ±1 or a = ±6. A Lax pair for (4) is given by [31] 

^ (5) 
M = — 4D^ — auDx — \aux I 

where D" denotes repeated application of Dj, (n times) and I is the identity oper- 
ator. Substituting C and M. into (2) yields 

Dli^ = (A - \au) i;, (6) 
DtV- = -403^; - au D^V - ^a^^xV'- (7) 
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The first equation is a Schrodinger equation for the eigenfunction ip with eigenvalue 
A and potential u{x,t). The second equation governs the time evolution of the 
eigenfunction. 
The compatibility condition for the above system is 

DtDl'^ - D^DiV = la {ut + auu^ + usx) ip = (8) 

where (6) and (7) arc used to eliminate Dttjj, D'^tp, D^tp and D^'^. Obviously, (6) 
and (7) will only be compatible on solutions of (4). 

The compatibility of (2) may be expressed directly in terms of the operators jC and 
M. Indeed 

Dt {jCi/j) = Ctip + >CDtV = ADtV 
where Ct'ip = Df {Cip) — CDftp. Using (2) we have 

AV' + ^M^p = XMip = MXtp = MC^P; 

that is, 

(A + CM- MC) ^ = 0. 

If this operator does not vanish then (2) is not involutive and so would have ad- 
ditional (nonlinear) constraints; that is, one cannot freely specify the initial data 
for (2). However, if this operator vanishes identically then Lax pair will not encode 
the original differential equation. Therefore, for a non-trivial Lax pair for (1), this 
operator must vanish only on solutions of (1). Hence 

Ct + [C,M] = 0, (9) 

where = denotes that this equality holds only on solutions to the original PDE (1). 
Here [C, M] = CM. — MC is the commutator of the operators and is the zero 
operator. Equation (9) is called the Lax equation. Note that Ct = [Dt,C] and so 
the Lax equation takes the form 

[Dt-M,C] = 0. (10) 
Example 2.2 Returning to (5) for the KdV equation, we have 

Ct = [Dt,C] = laut I, 
CM = - 4 - a(|u + ^u^ + (gaw^ + 2uxx) + {^auux + ^U3x) l) , 
MC = - 4 - a(^luDl + \ux + {\au^ + lu^x) + {\oiuUx + |u3c,) l) . 

Therefore, 

Ct + \C, M\ = \a {ut + auux + I, 
which is equivalent to (8). 



October 5, 2011 0:35 Applicable Analysis Hickmanetal-Applicable-Analysis-2011 

Applicable Analysis 5 

Various alternatives for the Lax pair operators exist. For example, one could define 
£ = C-X\ andM =_^M - Dt then (2) becomes Ctl; = and Mifj = 0. The Lax 
equation (9) is then [C, Ai] = 0. 

The order of the operator Ai may also be reduced. If the order of £ is £ then any 
terms in Ai may be rewritten in terms of derivatives of order at most £ — 1. 
Note that this occurs at the expense that the reduced operator M. has explicit 
dependency on the eigenvalue A. 

Example 2.3 For the KdV equation, using (5), we note that 

= [(A - lau) I] = (A - ^au) - ^au^ I (11) 
and so the third order operator M. is equivalent to 

;W = - (4A + lau) + lau^ I, (12) 

which is of first order but depends on A. 

There is a gauge freedom in the choice of a Lax pair. Suppose S be an arbitrary 
but invertible operator and let ^/j = Sip. Note that A^ = \Sil) = SjOtp = SCS~^'ip 
and 

However, [Dt,5]<S-^ = - 50*5"^ and so 50*5" V = SMS''^'^). Therefore, 

t = SCS-^, M=SMS-\ Dt = SDtS-^ (13) 
wiU satisfy [Dj -M,jC] = 0. 



3. Leix Pairs in Matrix Form 

In [2], Ablowitz et al. introduced a matrix formalism for Lax pairs. Their con- 
struction avoids the need to consider higher order Lax operators. They associated 
matrices X and T to the operators jC. and M respectively and considered the system 

14 

Dt^ = 

for an auxiliary vector function ^. Both X and T will be dependent on A. The 
number of components of ^ is determined by the order of £. When C is of order 
2, ^ has two components and X and T are 2x2 matrices. The compatibility 
condition for (14) is 

[Dt, D^l^- = Dt (X^-) - (T*) = (DtX) ^ + XD^^- - (D^T) ^ - TD^^ = 0; 

that is. 



(DtX-D^T+[X,T])* = 
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where [X, T] = XT — TX is the matrix commutator. Since we wish to freely specify 
initial data, the matrix Lax equation follows: 

DtX-D^T + [X,T] = 0. (15) 

Geometrically, we may interpret [44] the matrices X and T as defining a connection 
on a vector bundle over the base space given by (x, t) . ^' is an element of the vector 
bundle and parallel transport of this element on the base space is given by (14); 
that is, parallel transport along the x-direction is given by (Dj, — X) ^' = and 
along the t-direction by (Dj - T) ^ = 0. Note that 

[Dt - T, D, - X] = D^X - D,T + [X, T]. 

Hence, (15) states that the parallel transport of ^ is independent of the path taken. 
For this reason (15) is also known as the zero-curvature equation. In particular, we 
may rewrite (15) as [Vt, V^] = where = — T and Vx = D^; — X are the 
components of the covariant derivative induced by the connection. 

Example 3.1 It is well-known that 

X=f, °i I) (16) 



and 



^aux — 4A — ^au 

— 4A^ + laXu + joCt^u^ + \oiUxx — \oiU, 



form a Lax pair for the KdV equation since 



(17) 



DtX-D,T+[X,T] = -ia(^ 




ut + auux + U2,x 



which evaluates to the zero matrix on solutions of (4). 

Given an operator Lax pair, {C,M.), finding a corresponding matrix Lax pair, 
(X, T), is a straightforward but lengthy computation. As noted above in (11), the 
operator C allows one to write higher order derivatives in x in terms of lower order 
derivatives. Let the order of C be L Without loss of generality, we may assume 
that the leading coefficient of is 1. Thus, 

£ = D^ + /,_iDi-i + --- + /ol. (18) 

Let 



/ ^ \ 

^x^ 
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and so 



Do:* 



( 




VA-/0 



1 








1 





-/l -h ■ 

The computation of T is more elaborate. We have 



\ 


1 



(19) 



Note that, using (2) and (18), 

Df > = - + • • • + h D. + (/o - A) l) V. (20) 

Therefore any derivatives of ^ in a; of order I or higher are removed by applying 
(20) recursively. The result is a vector which will depend linearly on tp and its first 
i-1 derivatives; that is Dt* = T^'. 
We will only show the conversion explicitly for the KdV equation. 

Example 3.2 For the KdV equation, (19) reduces to (16). Using (12), we note 
that 

DfV' = \auxip - (4A + \au) D^^V, 
DtDj-V' = \oiUxx V' - Da;V - (4A + \au) Dlip. 

However, D^^ = (A — gon) 'ijj and so 

DtDj;^' = (-4A^ + \a\u + ^a^u^ + \auxx) - \oiUx Dxtp. 

Therefore, 



-aur 



-4A + -^aXu + j^a u + ^auxx —q 



which yields T in (17). 

There is gauge freedom in the construction of the matrix representation of a Lax 
pair. Let G be a non-singular matrix and * = G*. Using (14), 

D^* = GDx^ + (G) ^' = (GX + (G)) ^' = (GX + (G)) G'^^- = X# 
where 



X = GXG-^ + D,; (G) G-^ 



(21) 
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Similarly, Djl' = with 

T = GTG-^ + Dt (G) (22) 

Example 3.3 A second Lax pair [1] for the KdV equation consists of the complex 
matrices 



and 



~ \ -4P + ^au Uk^ — ^iaku + ^aux ' ' 

Based on (21) and (22), both Lax pairs are gauge equivalent with 

( -ik I 


where \ = —k"^. 



4. Dilation Invariance of Nonlinear Evolution Equations 

Crucial to the computation of Lax pairs is that (9) must hold on the PDE (1). This 
has the important consequence that any symmetry of the PDE, in particular the 
dilation symmetry, might be useful to find (9). In other words, we make an ansatz 
that the operators C and M inherit the scaling symmetry from the PDE. This, as 
it turns out, greatly simplifies the equations that determine the Lax pair. 

Example 4.1 The KdV equation (4) is dilation invariant under the scaling sym- 
metry 

{t,x,u) ^ {KT^t^KT^x^K^u), (25) 

where k is an arbitrary parameter. Indeed, after a change of variables with t = 
K~^t,x = K~^x,u = K^u and cancellation of a common factor the KdV for 
u{x, i) arises. This dilation symmetry can be expressed as tt ~ 9^ and dt ~ 
which means that u scales as two x— derivatives and the t-derivative scales as three 
x— derivatives. 

We define the weight, W, of a variable as the exponent of k in the scaling sym- 
metry. It is clear that, from (3), D^- ~ dx and Df ~ dt and so we extend weights 
to operators. For (25), W{x) = -1 (or Wid^) = l),W{t) = -3 (or W{dt) = 3) 
and W{u) = 2. However we could have scaled each of these weights (effectively 
K — )■ K*); that is, a scaling symmetry does not uniquely determine a set of weights. 
We use this freedom to set W{dx) = 1. 

The weight of a monomial is defined as the total weight of the monomial and 
will also be denoted by W. Such monomials may involve the independent and 
dependent variables and the operators dx, Dx-, dt and D^. In particular, note that 
W{Dx) = W{dx) and W{Dt) = W{dt). An expression (or equation) is uniform in 
weight if its monomial terms have equal weights. For example, (4) is uniform in 
weight since each of the three terms has weight 5. 
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The importance of uniformity of weight is demonstrated by the KdV equation 
(4). Prom (5), it is clear that C has weight 2 and M. has weight 3 since W{\) = 0. 
Therefore this Lax pair inherits the scaling symmetry. Also note that we must 
assign VF(A) = 2. The elements of the matrices X and T in (16) and (17) are also 
uniform in weight, albeit of different weights. This is also true for the matrices in 
(23) and (24). 



5. An algorithm for Computing Lax Pairs in Operator Form 

First, one computes the dilation symmetry of the nonlinear PDE and assigns 
weights such that W{dx) = 1. Next, one selects the weight £ of the operator jC. Since 
we have chosen W{dx) = 1, ( will also be the order of C The minimal weight for £. 
is the maximum weight of the dependent variables since wc want the Lax equation 
to depend non-trivially on the PDE. It is only through Ct that the t-derivatives 
of the dependent variables appear. Since Ct = [Dt,C], W{Ct) = W{C) + W{dt) 
and, from (9), W{Ct) = W{C) + W{M). Thus W{M) = W{dt). Since we require 
that to be a differential operator, we must choose a system of weights such that 
W{dt) is a non-negative integer. 
Thus the candidate Lax pair has the form 

C = Di + fe_,Di-' + ... + fo\, 
M = CmD^ + gm-i D--i + ---+5ol 

where m = W{dt). Now 

[C,M] = iiD,Cm)D^+^-^ + -- - . (26) 

Recall that we require that (9) only holds on solutions of the PDE. Thus the 
operator D^+^~^ could be reduced to an operator that depends on derivatives of 
order strictly less than i. However this would introduce the eigenvalue A and, in 
this case, there would be a term £X"^^^DxCm- This would be the only term in A™"^. 
Under the assumption that the operator C has a complete (and therefore infinite) 
set of eigenvalues, we conclude DxCm = 0; that IS J CfYi IS cl constant. Therefore (26) 
becomes 

[CM] = iiD,g^_, - mc^D,/,_i) D^+^-^ + ■■■ . 

After reduction, the coefficient of A"*"^ must vanish and so 

^D^gm-l - mCrrjL)xh-l = 0. (27) 

In this way, we obtain ra—1 equations for the unknown coefficient functions /«, gj 
that do not depend on the PDE. We call these equations kinematic constraints. 

The remaining (. components of (9) may involve terms from Ct- When t- 
derivatives do appear, we must use the PDE to remove them. We call these equa- 
tions dynamAcal equations. 

The kinematic constraints are easily solved. All have the form (27) and so may 
be solved in terms of the leading gj. Specifically, (27) yields 



9m— 1 — ^ (^mfl—\- 
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In this way, gj, j = 1, . . . , m — l, may be written in terms of go and fi. Note that 
this process does not depend on the PDE and so this rewrite is canonical. 

After the cUmination of the gj, the dynamical equations arc a set of ODEs for 
the remaining coefficient functions go and fi. We now make the ansatz that the 
terms in the operators C and M. have uniform weight; that is, fi has weight i — i 
and go has weight m. Thus, for example, go is assumed to be a linear combination 
(with unknown coefficients) of monomials of weight m. If, in addition, the weights 
of the dependent variables are positive and the monomials do not depend explicitly 
on X and t, the dynamical equations reduce to an finite system of algebraic equa- 
tions for the unknown coefficients. This algebraic system (which will also have any 
parameters that are present in the PDE) is then solved by Grobner basis methods. 

We now illustrate the steps of the algorithm for the KdV equation. 

5.1. Step 1: Computing the scaling symmetry 

The dilation symmetry of (4) can be readily computed by the requirement that (4) 
is uniform in weight. Indeed, setting W{dx) = 1, we have 



which yields W{u) = 2 and W{dt) = 3. This confirms (25). So, the requirement of 
uniformity in weight of a PDE allows one to compute the weights of the variables 
(and thus the scaling symmetry) with linear algebra. 

Dilation symmetries, which are special Lie-point symmetries, are common to 
many nonlinear PDEs. Needless to say, not every PDE is dilation invariant. How- 
ever non-uniform PDEs can be made uniform by extending the set of dependent 
variables with auxiliary parameters with appropriate weights. Upon completion of 
the computations one can set these parameters to 1. 

5.2. Step 2: Building a candidate Lax pair 

Since W{u) = 2, the minimal weight for C is 2. We therefore choose £ = 2. Note, 
that there is no monomial of weight 1 (any monomial that involves u must have 
weight at least 2). Therefore we could set f\ = 0. However the kinematic constraints 
and dynamical equations only depend on the choice of weight for C and W{dt). So, 
by not setting /i = at this stage, we will derive the dynamical equations for all 
Lax pairs with C of second order and W{dt) = 3. Our candidate Lax pair is 



The kinematic constraints (that is, the coefficients of D^ and in [£, A^]) give 



W{u) + W{dt) = 2W{u) + 1 = W{u) + 3, 



£ = Dl + fiDx + fo I, 
M = C3Dl+g2Dl+gi D^ + gol- 



(28) 



92 = fcs/i. 



gi = |c3 (2D^/i + /2 + 4/o) . 



Therefore, (28) must have the form 



£ = Dl + fiDx + fo I, 
M = C3Dl + |c3/i Dl + |c3 (2D^/i + /2 + 4/o) + go I. 



(29) 



October 5, 2011 0:35 Applicable Analysis Hickmanetal-Applicable-Analysis-2011 



Applicable Analysis 11 

The dynamical equations are 

Dt/i + 20^50 - |c3D^ (2D2/i + 12D^/o - ff + 12/i/o) = 0, 
Dtfo + Dlgo + fiD^go (30) 
- C3 (D3/o + IfiDlfo + |D,/iD,/o + IffDJo + IfoO^cfo) = 0. 

Equations (30) are the master equations for (29) . In other words, for any equation 
with a scahng symmetry such that W(dt) = 3, a Lax pair with a second order 
eigenvalue problem must have the form (29) where /o, /i, go and C3 satisfy (30). 
For the case of the KdV equation, we have /i = and so (30) reduces to 

2D,5o - fcsD^/o = 0, 
Otfo + Dlgo - C3 {Dlfo + |/oD^/o) = 0. 

Note that, in this case, the first dynamical equation does not depend on the PDE 
and so may be regarded as a kinematic constraint. This also implies that we cannot 
set C3 = in order to satisfy the kinematic constraint since this would result in a 
trivial dynamical equation Dj/o = 0. Therefore we have go = |c3Dj;/o and 

Dt/o - IcsDlfo - ic3/oD^/o = 0. (31) 
Since fo has weight 2 we must have 

fo = bou (32) 

with 60 7^ 0- 



5.3. Step 3: Computing the undetermined coefficients 

Equation (32) is substituted into the dynamical equation (31). We obtain 

hout - C3boU3x - ^csbluux = 0. 
Use (4) to replace ut by —{auux + u^x) and so 

60 [aUUx + U3x + {csU^x + ^CsboUUx) = 0. 

Equating to the coefficients of the monomials usx and uux, we have 

^'0 (1 + 2C3) = 0, 60 (a + §0360) = 0. 

Since 60 / 0, we must have C3 = —4 and 60 = Finally substitute fo = ^au, /i = 
0, go = —\aux and C3 = —4 into (29) to obtain the Lax pair (5). In this case the 
algebraic system was sufficiently simple that it did not require the machinery of 
Grobner bases to solve. 



6. Examples 



The method described in Section 5 can be applied to many nonlinear PDEs in 
(1 + 1) dimensions with polynomial nonlinear ities. In this section we apply the 
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method to scalar equations as well as systems. The featured scalar examples are 
the modified KdV and Boussinesq equations. The method is also illustrated for 
systems, including a coupled system of KdV equations due to Hirota & Satsuma 
and the Drinfel'd-Sokolov- Wilson system. 



6.1. The modified Korteweg-de Vries equation 

For the modified Korteweg-de Vries (mKdV) equation [1] for u{x,t), 

ut + au^Ux + = 0, (33) 

the weights may be chosen to be W{u) = 1, W{dx) = 1, and W{dt) = 3. Therefore 
the minimal weight for C is 1 and the weight of is 3. 
Suppose £. has order 1; that is, 



(34) 



jC = Dx + fo\, 
M = C3Dl+g2Dl+gi D^ + ^qI- 

The kinematic constraints are 

Dx92 - 3c3 D^/o = 0, D^gi - 3c3D^/o - 2g2Dxfo = 

and so 

92 = 3C3/0, gi = 3c3{Dxfo + /o)- 

The dynamical equation is 

Dtfo + D^go - csD^ (D^/q - 3/oD^/o - fi) = 0. 

In this case, we can replace go by go — c^ (D^/q — 3/oD2;/o — /q) to remove C3. Thus 
we can set C3 = without loss of generality. Therefore the Lax pair (34) is 

>C = D^ + /ol, 

^^^^ 

M=go\ 

with the dynamical (that is, master) equation 

DJo + 0^90 = 0. 
Now assume that /o has weight 1 and go has weight 3; that is 
fo = bou, go = aiu^ + a2UUx + a^u^x 
with 60 / 0. The dynamical equation becomes 

- 60 {au^Ux + uzx) + ^aiu^Ux + 02 {u^. + uuxx) + aa'^Sx = 0- 
Setting to zero the coefficients of the monomials usx, uUxx-> and u^Ux yields 
0-3 — bo = 0, a2 = 0, 3ai — boa = 0. 
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Since bo is free, we can set bo = 1 and (35) becomes 

'Cl = Da; + ul, 

(36) 

Ml = {^au^ + Ux^) I. 

The Lax equation for (36) is the conservation law form of (33). 

If we assume that C has weight 2, then, since W{dt) = 3, the Lax pair is given 
by (29) and the dynamical equations are given by (30) . In this case, a gauge on go 
will not remove C3 from the dynamical equations. The ansatz of uniform weights 
determines 

/i = bou, fo = biu^ + b2Ux, go = aiu^ + a2UUx + asu^x- 

The solution to resultant algebraic system has two branches; 

Case I: £ = DI + 2uDx + [u^ + u^) \, 

M = (lau^ + U:,:,) I. 
Case II: £ = + 2eu + i (6e\^ + 6eux + au^ ± V-6au^) I, 

M= -4:Dl- 12eu - (l2e\^ + 12eii^ + au'^ ± V-Qau^) 
— (4e^'u^ + |eaii^ + I2e^ uux ± ey/—6auux + S^Uxx 
+ auux ± ^V-6aUxx) I- 

The second case is a 1-parameter family (with parameter e) of Lax pairs. This 
family has the property — > — 4D^ for solutions of (33) such that n -> as 
|x| —7- 00. Wadati [38, 39] used this property to construct soliton solutions for (33) 
via the inverse scattering transform from Case II with e = 0. 
There is one Lax pair with C of order 3 which is 

jC = Dl + 3uDl + 3 {u^ + Ux) Dx + (u^ + 3uux + Uxx) I, 
M = {lau^ + Uxx) I- 

However 

jr = {Dx + u\f = Cl M = Mi 

where £1 and Aii are given by (36). 



6.2. The Boussinesq equation 
The wave equation, 

utt - Uxx + + Suuxx + (xu4^x = 0, (37) 

for u{x, t) with real parameter a, was proposed by Boussinesq [5] to describe surface 
waves in shallow water [1] . This equation may be rewritten as a system 



ut-Vx = 0, vt-Ux + 3uux + au^x = 0, 



(38) 
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where v{x,t) is an auxiliary dependent variable. We could proceed following the 
above examples modulo the discussion on scaling below. However we will work 
directly with (37). First note that (37) can be written as a conservation law, 

Dtut = Da; (ux - Suux - au3x) ■ (39) 

However it is not uniform in weight (the monomials Ux and uux cannot have the 
same weight unless W{u) = 0). This problem may be circumvented by the intro- 
duction of an auxiliary parameter j3 with an associated weight. Thus, we replace 
(39) by 

DtUt = Dx {I3ux - 3uux - ausx) (40) 

and assign weights 

W{dx) = 1, W{dt) = 2, W{u) = 2 and W{/3) = 2. 

This example demonstrates that a PDE, which is not dilation invariant, may be 
expressed in a dilation invariant form by the introduction of auxiliary parameters 
with weights. Upon completion of the computations, one can set each of these 
parameters to 1. 

The "fundamental" dependent variable in this equation is ut (that is, iiu = ^tUt)- 
However the Boussinesq equation is a conservation law. Therefore we can use D~^ut 
as a fundamental variable; in other words, rewrite (40) as 

Dt [Dx^ut) = /3ux - Suux - ausx- 

Therefore the minimal order of C is 3. Since there are no weight 1 monomials, the 
minimal order Lax pair is given by 

ji: = Dl + hDx + M, 

M = C2Dl + go\. 
The kinematic constraint is SOxgo = 2c2D-r/i which gives 

90 = ic2/i + Co/3. 

The extra term arises since the constant /3 has weight 2. The dynamical equations 
become 

Dth = C2 (2D,/o - Dlh) , (41) 
Di/o = C2 (D^/o - plf, - IhDxh) . 

Since it is uu and not Ut that appears in (40), derivatives of u that have at most 
only one t-derivative may be freely specified in the initial data for (40). Thus 
Uf, Uxt, Uxxt, ■ ■ ■ may appear in monomials used to construct /o and /i. In addition, 
/o only appears in the right hand side of the dynamical equations in the form Dxfo- 
Thus a non-local term of the form D~^f, where / has weight 4,may appear in /q. 
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Consequently, the weight ansatz imphes 
/i = aiu + 02/3, 

/o = CLsUx + D^:^ [a^u^ + a^Pu + UQUt + 07^^) . 
Under this ansatz, the dynamical equation (41) becomes the kinematic constraint 

aiut = 2c2 {a^Uxx + a^u^ + a^^u + aeUt + ay/?^ - ^aiUxx) ■ (42) 

For a non-trivial solution, we require ai 7^ and so C2 7^ 0. Thus 04 = as = 07 = 0, 
ai = 2c2aQ and 03 = C2aQ. Hence, 

Dt/o = asUxt + aeD~^utt = ae ic2Uxt + Pux - 3uux - au3x) 
on solutions of (40) . Therefore the remaining dynamical equation is 

ae {c2Uxt + PUx - ^UUx - "'"3a;) = C2a6{uxt - - ^C2Ux{2c2aQU + 02/3)). 

The solution of this equation gives the Lax pair 

>C = D3 + i- (3u - /3) + ^ {aux ± D-^^t) I, 

2a 

The subcase /3 = 1 is a Lax pair for (37) [43]. Note that the auxiliary variable v 
that was introduced to obtain the system (38) is, in fact, D~^ut. A Lax pair for 
the system is obtained by replacing ^x^ut by v. 

6.3. The coupled Korteweg-de Vries equations 
The coupled Korteweg-de Vries (cKdV) equations [21], 

ut - 6/3uUx + 6vvx - Pusx = 0, vt + Suvx + V3x = (43) 

where /3 is a nonzero parameter, describes interactions of two waves with different 
dispersion relations. System (43) is known in the literature as the Hirota-Satsuma 
system. It is completely integrable [1, 21] when /3 = ^. 

We assign the weights W{dx) = 1, W{dt) = 3, W{u) = 2 and W{v) = 2. Since 
W{dt) = 3 and there are no weight 1 monomials, the dynamical equations are 
precisely those of the KdV equation (4). The minimal weight for C is 2 and so the 
dynamical equation is (31) with 

/o = bou + biv. 

The resultant equations have only a trivial solution. Likewise, there are also no 
Lax pairs with £, of third order. 
For the case W{jC) = 4, the candidate Lax pair, after the kinematic constraints 
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have been solved, is 

C = Dt + f2Dl + hD, + fo\, 
M = C3Dl + f C3/2 + |c3 (2/1 - D^/2) I 

with dynamical equations 

Dt/2 + \c3 {GDlfi - Dlf2 + 3/20^/2 - 12D^/o) = 0, 
Dt/i + (8D^/i - 3D^/2 + 3/2D^/i + 3/iD^/2 - 12D2/o) = 0, 
Dt/o + |c3 (6D^/i - 3D^/2 + 6/2D^/i - 3/20^/2 

+ 6/iD^/i - 3/iD2/i - 8D3/o - 6/2D^/o) = 0. 

A non-trivial solution exists only when P = ^.In this case the solution is 

£ = D^^ + 2uDl + 2{u^- v^) + [u^ - v'^ + - V2x) I, 
X = 2 + 3n + I (uj; - 2vx) I. 

Note that this Lax pair is given in [9, 42]. 



6.4. The Drinfel'd-Sokolov- Wilson equations 

We consider a one-parameter family of the Drinfel'd-Sokolov- Wilson (DSW) equa- 
tions 

ut + Svvx = 0, vt + 2uvx + auxv + 2vzx = 0, (44) 

where a is a nonzero parameter. The system with a = 1 was first proposed by 
Drinfel'd and Sokolov [12, 13] and Wilson [42]. It can be obtained [26] as a reduction 
of the Kadomtsev-Petviashvili equation (a two-dimensional version of the KdV 
equation). 

We may assign weights W(dt) = 3, W{u) = 2 and W{v) = 2. Thus we have the 
same weight system as the KdV equation and the coupled system above. There are 
no Lax pairs for this system with jC of order 2, 3, 4 or 5 which inherit this scaling 
symmetry. 

There are no solutions for W{£,) = 6 except for a = 1. In that case 
£ = Dl + 2uDl + {4ux - 3vx) + (| {u2x - V2x) + f ') 

+ (i (uSx - V3x) + 2 {uUx - VVx) + UxV - UVx) Dj, 

+ i (2 (U^x - VAx) + 2{u + v) {u2x - V2x) + - vl) I, 

M = Dl + uDx-\{'ivx-Ux)\. 

The DSW system (44) is known to have infinitely many conservation laws when 
a = 1 and is completely integrable [12, 42]. This is the Lax pair given in [42]. 
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Parameter 

ratios 

a f3 
^2' ^) 


Commonly 
used values 

(«,A7) 


Equation name 


Lax pair 
references 




(30, 20, 10), (120, 40, 20), 
(270,60,30) 


Lax 


[31] 


(ii) 


(5,5,5), (180,30,30), 
(45, 15, 15) 


Sawada-Kotera 


[10, 15] 
[27, 36] 


/l 5\ 
(5' 2) 


(zU, zo, lUj 


Kaup-Kupershmidt 


[15, 27] 


(1,2) 


(2, 6, 3), (72, 36, 18) 


Ito 





Special cases of the fifth-order family (45) of KdV-type equations. 



7. Classification of fifth-order KdV-like equations 

For an equation or system with parameters, the algorithm described in Section 5 
can be used to find the necessary conditions on the parameters so that the PDE 
has a Lax pair and therefore may be completely integrable. 
Consider the family of fifth-order KdV-type equations with three parameters 

Ut + aU^Ux + Pu^Uxx + 7UU3X + Ubx = 0. (45) 

This family includes several well-known completely integrable equations. 

Replacing by ^, it is clear that only the ratios ^ and ^ matter. We continue 
with (45) for an easier comparison with the results in the literature. The special 
cases of (45) shown in the table are extensively discussed in the literature [18, 22, 
29, 36]. Only the first three equations are completely integrable. Ito's equation is 
not completely integrable but has an unusual set of conservation laws [18, 23]. 

We wish to determine which members, if any, have a Lax pair, with C of second 
order, which may then be amenable to a (standard) inverse scattering transform 
based on a quadratic eigenvalue problem. Kaup [27] has considered this family with 
respect to a cubic eigenvalue problem (that is £ is a third order operator). 
Uniformity of weights implies that W(dt) = 5. The first kinematic constraint is 

2Da;C/4 - bcsDxfi = 0. 

In this situation we may need or wish to give the parameters weight. If this is the 
case, there may be a weight 1 constant. Therefore the solution to this constraint is 

94 = ic5/i + C4 

where C4 is a constant of integration of weight 1. After all four kinematic constraints 
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are solved, we have 

C = Dl + hD^ + fo I, 

M = C5Dl + + C4) + (|c5(6D,/i + 3/2 + 4/o) + 2C4/1 + C3) 

+ (^C5(12D2/i + + 12D^/o + ff + 12/i/o) + C4(2D^/i + + 2/o) 

+ fcs/i + C2) D2 + (^C5(48/oD,/i + 96/iD,/o + 4D^ff + 24f!fo 

+ 20(D,/i)2 + 40/iD2/i + 24D3/i + SOD^/o + 48/^ - ff) 

+ C4(2D^/o + /iD^/i + Dlfi + 2/1/0) + |c3(4/o + /f + 2D^/i) 

+ C2/1 +cij Dj; + 5-ol- 

For brevity, we omit the two dynamical equations. 

We now need to make a choice for the remaining weights. The obvious choice 
is W(u) = 2. With this choice, a, (3 and 7 all have weight zero. Therefore the 
constants of integration, ci, . . . ,04, are all zero which simplifies the Lax pair dra- 
matically. The only non-trivial solution occurs when a = ^7^, /3 = 27 which 
corresponds to the fifth-order KdV equation due to Lax [31]. The Lax pair is 

£ = D2 + ^7«|, 

M= - 16 - A-fu Dl - 67^x^ - 7 (5?x^^ + ^-iu^) D^^ - 7 (|w3^ + ^-tuu:,) I. 

In [7], it is shown that all scalar evolution equations, which have the form of a 
conservation law, have a Lax pair with a second order C. Such Lax pairs can be 
found if we choose the weights W{u) = 1, W{a) = 2, W{I3) = 1 and ^(7) = 1. 
In this case, the constants of integration may not be zero. 

For weight 3, we obtain non-trivial solutions for the Sawada-Kotera and Kaup- 
Kupershmidt equations. For the Sawada-Kotera equation (a = ^7^, 
beta = 7) [8, 37] there are two solutions: 

Case I: £. = 01 + ^ju D^, 

M = 9Dl + 3-fuDl + 3ju^ 0^+7 [l^u'^ + 2n2^) 
Case II: £ = 01 + ^7^ + ^jUx X, 

M = 9Dl + 3-fuDl + 6jux 0^+7 (17^^ + 5u2x) 

+ 27 i^lUUx + USx) I- 

The first Lax pair is given in [10, 15]. 
The Kaup-Kupershmidt equation (a = ^7^, P = I7) [27] has a single Lax pair 

C = DI + Iju D^ + ^7u^X, 

which is given in [15]. We cannot find any reference to the second Lax pair in the 
literature. 
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8. Conclusions 

The algorithm described in this paper can be easily implemented in Mathematica 
or Maple. Therein lies its value; it gives a quick method to "test the waters" with 
a new system of equations. If the equations do not have a scaling symmetry, or 
the scaling symmetry does not yield a Lax pair, then weighted parameters may be 
introduced. A new assignment of weights may lead to a Lax pair. The "focus on the 
equation" of our approach is a marked contrast to existing work where integrable 
systems have, by and large, been developed from geometric considerations (for 
example, the use of the Drinfel'd-Sokolov construction [12, 13]). 

The question of an appropriate system of weights is non-trivial. As noted above, 
all members of the three parameter family of fifth-order KdV-like equations have 
a Lax pair with a second order C. The question naturally arises as to which of 
these Lax pairs are useful. Can, for example, a criteria be developed that would 
guarantee an infinite family of local conservations laws? 
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